Fast Neutron Detection with ^Li-loaded Liquid Scintillator 

B. M. Fisher^'\ J. N. Abdurashitov^, K. J. Coakley^ V. N. Gavrin^, D. M. Gilliam^, J. S. Nico^, A. A. Shikhin^, 

A. K. Thompson^, D. F. Vecchia^ V. E. Yants^ 

^National Institute of Standards and Technology, Gaithershurg, MD 20899 USA 
^Institute for Nuclear Research, Russian Academy of Sciences, Moscow, 117312 Russia 
^National Institute of Standards and Technology, Boulder, CO 80305 USA 



Abstract 



o 

(N 

(N 
(N 



• 1—1 

O 

• 1—1 



We report on the development of a fast neutron detector using a liquid scintillator doped with enriched ^Li. The lithium 
was introduced in the form of an aqueous LiCl micro-emulsion with a di-isopropylnaphthalene-based liquid scintillator. A 
^Li concentration of 0.15 % by weight was obtained. A 125 mL glass cell was filled with the scintillator and irradiated with 
fission-source neutrons. Fast neutrons may produce recoil protons in the scintillator, and those neutrons that thermalize 
within the detector volume can be captured on the ^Li. The energy of the neutron may be determined by the light output 
from recoiling protons, and the capture of the delayed thermal neutron reduces background events. In this paper, we 
discuss the development of this ^Li-loaded liquid scintillator, demonstrate the operation of it in a detector, and compare 
its efficiency and capture lifetime with Monte Carlo simulations. Data from a boron-loaded plastic scintillator were 
acquired for comparison. We also present a pulse-shape discrimination method for differentiating between electronic and 
nuclear recoil events based on the Matusita distance between a normalized observed waveform and nuclear and electronic 
recoil template waveforms. The details of the measurements are discussed along with specifics of the data analysis and 
its comparison with the Monte Carlo simulation. 
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1. Overview 

Fast neutrons may be produced through several mecha- 
nisms. Naturally occurring isotopes in the ^^^U and ^32^^]^ 
decay chains generate fast neutrons via spontaneous fis- 
sion, and (a,n) reactions create neutrons in the 1 MeV to 
15 MeV range. Particle accelerators produce higher en- 
ergy fast neutrons through a variety of reactions, and neu- 
trons of very high energies are generated from cosmic-ray 
muon-induced spallation reactions [1 . Improvements in 
the ability to detect and characterize these neutrons are 
of interest to diverse research interests such as fundamen- 
tal physics, neutron dosimetry, and detection of low-level 
neutron emissions. 

In the area of basic physics research, many classes of nu- 
clear, particle, and astrophysics experiments must be per- 
formed in underground laboratories to reduce the back- 
grounds generated by naturally occurring radioactivity 
and cosmic rays [2 . Although experimenters operate their 
detectors deep underground and go to great lengths to op- 
timize radiation shielding, the experiments have become 
so sensitive that even these efforts may not be sufficient. 
The characterization of the fast neutron ffuence has be- 
come a critical issue for experiments that require these ex- 



treme low-background environments, such as neutrinoless 
double-beta decay [SlIHIS], dark matter searches O [3 [3 [9] , 
and solar neutrino experiments [10, 11, 12^ ITSj [HI [15]. In 
some experiments, fast neutrons may be the dominant and 
potentially irreducible background, thus necessitating pre- 
cise information about the fast neutron ffuence and energy 
spectrum. The most reasonable approach to addressing 
the problem is through the complete characterization of 
the neutrons through both site-speciffc measurement [16] 
and benchmarking of simulation codes [171 HI] • 
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The health physics community is another area where 
improved fast neutron detection and spectroscopy are 
needed [19 . Existing spectrometers allow the determina- 
tion of ffuence spectra (and dose) for low energy neutron 
sources such as ^^^Cf and D2O moderated ^^^Cf, but be- 
gin to have difficulty with higher energy sources such as 
^^^Am-Be(a, n) or ^^^Am-B(a, n). Current spectrometers 
fail almost completely for determining neutron ffelds with 
energies of tens of MeV and may require multiple measure- 
ments with different detectors and complicated unfolding 
procedures. This need has only grown due to the increased 
use of 14 MeV neutron generators in interdiction and in- 
spection technologies. Without good knowledge of these 
neutron spectra, health physicists cannot accurately cal- 
ibrate radiation protection instruments and dosimeters, 
which may result in incorrect determination of dose re- 
ceived by radiation workers. 



Preprint submitted to Nuclear Instruments and Methods A 



June 23, 2011 



An improved fast neutron detector has direct applica- 
tion to the detection of low fluence rates of fast neutrons, 
such as from fissile material. The technological challenges 
and requirements for measuring the fast neutron fluence 
in the underground environment are very similar to this 
problem. A highly efficient neutron spectrometer with rea- 
sonable resolution would be capable of detecting low-level 
neutron signals from a number of sources. 

In this paper, we present the developmental work per- 
formed using a small detector filled with an enriched ^Li 
liquid scintillator. General methods of neutron detection 
and spectroscopy may be found elsewhere [20l [21]. Sec- 
tion |2] discusses capture-gated neutron detection and the 
fabrication of the ^Li-loaded liquid scintillator. In Sec- 
tion [3] the results of measurements with a test detector of 
^ Li- loaded liquid scintillator and a B-loaded plastic scintil- 
lator are presented along the analysis methods. Section [4] 
presents a pulse-shape discrimination method for differen- 
tiating between electronic and nuclear recoil events. 

2. Neutron Detection with ^Li Scintillator 

2.1. Capture- Gated Neutron Detection 

The detection method used in this work is known as 
capture-gated neutron spectroscopy [22l |23l [24l |25]. As 
illustrated in Figure [l] an incident fast neutron prefer- 
entially scatters from a proton in a scintillation medium. 
The proton recoils with approximately half (on average) 
of the neutron's energy, producing scintillation light. The 
neutron then scatters off another proton, producing an- 
other proton recoil, and continues to scatter until it leaves 
the detector volume or is captured, typically at thermal 
energies. A large fraction of the neutrons that lose all 
their energy within the scintillator volume and then cap- 
ture will have lost most of their energy through neutron- 
proton scattering, and the majority of the detected light 
will come from those neutron-proton events. There are 
other mechanisms that may generate light, in particular 
scattering from carbon [26 , but these are small effects for 
the measurements performed here. 

An organic scintillator provides a high density of pro- 
tons that will efficiently moderate the neutron; a neutron 
with an energy of a few MeV loses 90 % of its energy in 
this manner during the first 10 ns in the scintillator. The 
thermalized neutron diffuses within the bulk of the scintil- 
lator and may capture on an isotope that was introduced 
because of its large neutron capture cross section and un- 
ambiguous signature of the neutron capture. The time to 
capture occurs within approximately 100 /is and depends 
upon upon the size and geometry of the detector and the 
concentration and neutron capture cross section of the iso- 
tope used as the dopant. The delayed-coincidence signa- 
ture indicates that the fast neutron gave up essentially all 
of its energy within the scintillator and provides a method 
of rejecting of uncorrelated backgrounds. 

One must consider several factors when selecting which 
capture isotope is appropriate for the application. Typical 




Figure 1: An illustration of the principle of capture-gated detec- 
tion. A fast neutron impinges on the detector. It rapidly gives up 
its energy through nuclear collisions, primarily with protons, in the 
moderation process. The thermalized neutron diffuses in the medium 
until it is captured on a material with a high capture cross section. 

elements are gadolinium, boron, and lithium. Gadolinium 
is commercially available and has a very large neutron cap- 
ture cross section, but it produces gamma rays of various 
energies (up to 9 MeV) that may be difficult to detect with 
small volume detectors. In addition, the Gd gamma rays 
cannot be distinguished from background gamma rays by 
means of pulse shape discrimination. With neutron cap- 
ture on ^^B, the ^Li can be formed in a ground or excited 
state, leading to two reaction branches: 

+ n ^ ^Li* (0.84 MeV) + a (1.47 MeV) (1) 
^Li* -^^ Li + 7 (0.477 MeV) 
^°B + n ^ '^Li (1.01 MeV) + a (1.78 MeV). (2) 

The branching ratio for the first reaction is 93.7% and 
6.3 % for the second. The cross section is also large, and 
the alpha particle makes efficient thermal detection pos- 
sible with a small- volume detector. Boron- loaded detec- 
tors are expensive, which tends to make scaling to larger- 
volume detectors cost-prohibitive. 
A good alternative to ^^B is ^Li 

^U^n^t (2.05 MeV) + a (2.73 MeV). (3) 

The use of enriched ^Li as the dopant has the advantage 
of a large Q-value and the production of two energetic 
charged-particles. The light yield of the 2-MeV triton is 
nearly a factor of 10 higher than that of the 1.5-MeV alpha 
from neutron capture on ^^B. The neutron capture signa- 
ture from ^Li is well-separated from the noise and most 
background sources. There is no concern about the en- 
ergy leaving the scintillator as there is with dopants that 
produce gamma rays, and pulse shape discrimination be- 
tween gammas and the capture products is possible. The 
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Isotope Scintillator cjn /(%) Volume H/C 



"^i 
B 



liquid 
plastic 



941b 
3842 b 



0.15 
1.1 



125 cm'^ 
103 cm^ 



1.5 
1.1 



Table 1: Properties of the two scintillator materials used in testing. 
(Jn is the 2200 m/s neutron capture cross section; / is the fractional 
weight of the isotope in the scintillator; and H/C is the hydrogen to 
carbon ratio of the scintillator. 
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capture cross section is high, but it less than for Gd and 

The capture-gated detection method using ^Li has been 
implemented in different ways. A 5-L detector was con- 
structed for measurements in the Gran Sasso Labora- 
tory [27^, and an 8-L detector was used to measure the 
fast neutron background at the Modane Underground Lab- 
oratory in Modane, France ^28^. Both detectors used the 
commercial ^Li-doped liquid scintillator NE320 [29]. In a 
different approach, a plastic scintillator impregnated with 
lithium gadolinium borate crystals was developed and was 
demonstrated to detect neutrons over a wide range of en- 
ergies [24] . Studies of the response of B-loaded liquid scin- 
tillator to monoenergetic neutrons have also been carried 
out [30]. 

2.2. ^Li-loaded Liquid Scintillator 

Natural lithium for use as a neutron capture material 
was introduced into organic liquid scintillators more than 
50 years ago [SU |^ . More recently a solution made using 
^Li-salicylate dissolved in a toluene-methanol solvent [33 
was developed. Scintillators with enriched ^Li were de- 
veloped later, and in the 1980s Nuclear Enterprises manu- 
factured a lithium-loaded pseudocumene-based scintillator 
(NE320 [34 ), which was used to make several measure- 
ments of neutron backgrounds in underground laboratories 
[29) [271 EH] . This scintillator is no longer being manufac- 
tured, and to the authors' knowledge there are currently 
no vendors of scintillator doped with enriched ^Li. Hence, 
detectors incorporating ^Li have not been made for many 
years due to the lack of commercial availability of a scin- 
tillator. 

Given the possible applications for a ^Li-doped scin- 
tillator, we have investigated approaches for introducing 
lithium compounds into liquid scintillators. After explor- 
ing methods for incorporating organolithium compounds 
directly into an organic scintillator, we investigated tech- 
nology used in liquid-scintillation counting in the life sci- 
ences. Radioactive samples used in biochemistry are often 
low-energy beta emitters in an aqueous solution, such as 
tritium in urine. The very short range of these low-energy 
/3-particles make traditional counting very difficult. The 
liquid scintillators used in this field have additional sur- 
factants that allow for the creation of micro-emulsions be- 
tween the aqueous solution and the organic liquid scintil- 
lator. This permits the counting sample to be uniformly 



Figure 2: Block diagram of the data acquisition electronics. The two 
inputs to the digitizer allow for different range settings. 

distributed throughout the scintillator, greatly increasing 
its counting efficiency for low energy /3-particles. 

A critical property of the desired organic scintilla- 
tor was its ability to distribute and suspend a lithium- 
containing compound uniformly throughout its volume. 
We found a scintillator manufactured by Zinsser Ana- 
lytic [34] composed of a solvent (very high-purity di- 
isopropylnaphthalene) with added surfactants allowing in- 
troduction of aqueous solutions up to 40 % water by vol- 
ume. Di-isopropylnaphthalene is considered a safer solvent 
than the more common pseudocumene-based liquid scin- 
tillators. It has a higher flashpoint (approximately 150 °C) 
and is biodegradable, allowing for safer handling and op- 
eration. 

A 10- molar aqueous solution of ^Li-enriched LiCl was 
created by reacting enriched Li2C03 (95 % ^Li by weight) 
with concentrated hydrochloric acid, boiling off the excess 
acid, and dissolving the remaining dried LiCl salt in pure 
de-ionized water. Mixing this LiCl (aqueous) solution with 
the liquid scintillator, we achieved ^Li concentrations of 
approximately 0.15% by weight [35]. This concentration 
is comparable to what had been commercially available 
and used previously [28, 27 . This scintillator mixture was 
found to have good light properties, and no precipitation 
of LiCl was seen during the duration of testing. 

3. Detector Characterization 

3.1. ^Li-loaded Liquid Scintillator Detector 

A 125-mL cell was filled with this scintillator and in- 
serted into fast neutron fields to demonstrate its effective- 
ness as a neutron detector. Some characteristics of the cell 
and scintillator are given in Table [1] A 5.08 -cm diameter 
photomultiplier tube (PMT, Burle 8850) was mounted on 
one face of the 5.08-cm diameter by 5.08-cm long cylin- 
drical cell, and the remaining walls were coated with a 
diffuse, white reflector. The logic of the data acquisition 
electronics is shown in the block diagram of Fig. |2j The 
PMT signal was amplified and split into two signals. One 
output was delayed 64 ns to avoid self-triggering, and it 
became the start signal in a time-to- amplitude converter. 
To further reduce the probability of self-triggering, the dis- 
criminator was vetoed for 500 ns after an event. The trig- 
ger was formed on any event above a threshold followed by 
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Figure 3: Plot of the digitized PMT pulse of the detector response to 
a fast neutron within the scintillator. A neutron enters the test cell 
and initially scatters off the protons in the scintillator. After ther- 
malizing within the scintillator, it captures on a ^Li, in this example 
approximately 1.4 /is later. 



a second event that occurred within 40 /is; the waveform 
containing both events was digitized and recorded to disk. 
The data acquisition system was a single 2-channel 2 GHz 
digital oscilloscope card with 8-bit resolution. 

Figure [3] gives a sample trace showing the initial pro- 
ton recoil of a fast-neutron off hydrogen in the scintillator 
followed by the ^Li-capture pulse. All of the analysis was 
performed on these digitized waveforms. Each waveform 
was analyzed to determine the energy of each pulse and the 
relative timing of the two events. The pulse offset was de- 
termined by averaging 200 ns of the baseline near the stop 
event. The pulse height corresponds to the maximum volt- 
age of the pulse, and the time of the pulse was obtained 
from the timing channel corresponding to the position of 
the maximum pulse height. Histograms of these param- 
eters from the measured data were constructed from the 
analysis of all events. 

The light yield of the scintillator is nonlinear as a func- 
tion of proton recoil energy. Thus, the sum of all the light 
from individual recoil events is not proportional to the in- 
cident neutron energy. This results in the degradation of 
the energy resolution. By segmenting the detector, one 
can ensure that the probability of having more than one 
energy deposit in any segment is very small and thus re- 
construct the initial neutron energy by summing the light 
output of the individual segments [36 . This produces a 
detection scheme that can still be efficient while achiev- 
ing good energy resolution and suppression of uncorrelated 
backgrounds. There is an effort in progress to construct a 
16-channel pilot spectrometer based on this idea [37] . 

The test cell was exposed to fission neutrons from ^^^Cf 
sources and 2.5-MeV monoenergetic neutrons from a com- 
mercial neutron generator. Both the sources and generator 
are maintained at the NIST Californium Neutron Irradia- 
tion Facility (CNIF) [38 . The top plot of Figure [I] shows 
the energy spectrum for both the start and stop pulses of 
the trigger when the detector was irradiated with neutrons 
from a ^^^Cf source. The start events correspond primarily 
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Figure 4: Results from a demonstration of the ^Li-loaded liquid scin- 
tillator cell irradiated by a ^^^Cf source. The top plot shows the 
energy distribution of the start events (primarily recoil events) and 
stop events (primarily capture on ^Li). The peak in the spectrum 
of stop events corresponds to the total energy of the ^Li reaction 
products. The bottom plot shows the delay time (between start and 
stop events) distribution for event sequences where the stop event 
produces light consistent with neutron capture on ^Li. The error 
bars represent plus and minus the combined standard uncertainty. 



to proton recoil events while the stop events primarily cor- 
respond to the capture of the thermalized neutron on ^Li. 
The peak in the capture spectrum is due to the reaction 
products from the ^Li(n,a)^H reaction; most of the scin- 
tillation light comes from the triton. The peak is clearly 
separated from the noise, and the signal-to-background is 
more than 10:1. The energy spectra were calibrated using 
the Compton edges of gamma sources such as ^^^Ba and 
^^^Cs, as there is little photopeak detection in the liquid 
scintillator. Note that the energy scale for both plots is 
given in the electron-equivalent energy. The timing data 
were fit to a two-component exponential with a constant 
background, as shown in the bottom plot of Figure [4j The 
error bars represent plus and minus the combined standard 
uncertainty. A long component arises from the capture of 
thermalized neutrons as they diffuse within the bulk of the 
scintillator. The simulation indicates that at times below 
1 /is there is a population of partially thermalized neutrons 
that can still capture on ^Li or ^^B. These events give rise 
to a second, faster component in the timing spectrum. 

Because the data were digitized event-by-event, recoil 
events that do not correspond to a valid neutron capture 
can be rejected in analysis. Two important cuts are the 
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energy of neutron capture and the timing of the capture. 
An energy cut can be made around the ^Li capture peak 
(shown in Fig. [4| to reject with high probabiUty recoil 
events that do not have a vahd capture associated with 
them. The capture peak was fit to a gaussian function, 
and the energy window was chosen to be ±3 sigma around 
the peak channel of the distribution. The bottom plot 
of Fig. [4] shows the time distribution of capture times for 
the thermalizing neutron after the energy cut was made 
to keep only events with neutron captures on ^Li. One 
can use the timing spectrum to perform another cut. At 
long times, the correlated capture events are gone, and 
one is left with the uncorrelated background. The spec- 
trum of these background events can be subtracted from 
energy spectrum of the recoil events. Typically, the signal 
window was 0.5 /is to 10 /is and the background window 
was 10 /is to 20 /is. The relative size of the windows could 
vary slightly from run to run depending on the signal-to- 
background in the particular time spectrum. 

Data were also acquired in a similar data manner us- 
ing plastic scintillator loaded with 5 % natural boron (BC- 
454), corresponding to about 1 % ^^B. Some its character- 
istics are given in Table [T] The scintillator was a 5.08 cm 
by 5.08 cm right cylinder coupled to a 5.08 cm diameter 
PMT (Burle 8850). The setup, data acquisition, and anal- 
ysis were carried out in the same manner as that done for 
the ^Li-loaded cell. These data were useful for compari- 
son with the ^Li data and also as another check on the 
computer simulations. 

The top plot of Figure |5] shows the energy spectrum for 
both the start and stop events when the detector was ir- 
radiated with neutrons from a ^^^Cf source. The energy 
scale was calibrated using a ^^^Ba source and the Comp- 
ton edge of the 477 keV photon, which is seen in the stop 
spectrum. The number of events for the gamma is less 
than for the alphas due to the lower efficiency for stopping 
the gammas in the small volume of plastic. The bottom 
plot of Figure |5] shows the time distribution of neutron 
captures on the ^^B. As with the ^Li cell, the distribu- 
tion was made after an energy cut on alpha peak in the 
stop spectrum, thus requiring that there was a valid neu- 
tron capture. The distribution of delayed-capture times 
fits well to the 2-component exponential function. 
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Figure 5: The top plot shows the energy distribution of the start and 
stop events for the irradiation of plastic scintillator by neutron from 
252 Q£ rpj^g large peak in the spectrum of stop events comes from the 
alpha particle and the smaller peak to the right is the Compton edge 
of the 477 keV gamma. The bottom plot shows the time distribution 
of all events that have a valid neutron capture on ^°B. 



the detected rate, detector composition, background sub- 
traction, and dead time. Table [2] gives some of the rates 
and run parameters associated with the data set. 

When irradiated with a source, the neutron rate Rn is 
the difference between the total measured rate Rtot and 
the uncorrelated background rate Rbkgd that is associated 
with the source. The numerical values are obtained after 



applying both the energy and timing cuts of Section |3.1 
The neutron rate can also be expressed as 



Rn — Rti 



R. 



bkgd 



eAcf^ + Rrr + Rb, 



(4) 



3.2. Efficiency Measurements 

Both the ^Li cell and the B-plastic scintillator were irra- 
diated by ^^^Cf sources of known activity in the CNIF, and 
measurements of their efficiencies for fast neutron detec- 
tion were performed. Irradiations were carried out using 
two source activities and two gamma shielding configura- 
tions, and the measurements were done at several source- 
detector distances. For the gamma shielding, there was 
either a thin (6 mm) annulus of lead around the detector 
or there was no shielding at all. The purpose in using 
several configurations was to test the robustness of the 
analysis and simulation under varying conditions, such as 



where e is the detector efficiency, Aci is the neutron source 
activity, and Q is the solid angle subtended by the detec- 
tor from the source. Rrr is the rate in the detector due 
to room-return neutrons (i.e., those source neutrons that 
scatter from the surrounding environment into the detec- 
tor), and is the background rate in the detector when 
there is no source present. The measured rates must be 
obtained above a low-energy threshold. The values were 
chosen to be safely above the hardware threshold and the 
detector noise, and they are listed for each configuration 
in Table [1 

The background rate R^ is negligible compared with the 
room return and can be ignored for these measurements. 



5 





Pb shield 


Activity 


Number of 


Range 


RtoI 


Rbkgd 


Threshold 






(Bq) 


Positions 


(cm) 


(s-1) 


(s-1) 


(keV) 


^Li liquid 


off 


1.2 X 10^ 


7 


107 - 230 


3.90 


3.77 


75 


scintillator 


on 


1.2 X 10^ 


6 


86 - 230 


3.72 


3.40 


76 




off 


1.0 X 10^ 


6 


8-33 


0.308 


0.258 


86 




on 


1.0 X 10^ 


6 


9-30 


0.0530 


0.023 


93 


B-plastic 


off 


1.2 X 10^ 


7 


81 - 231 


5.07 


4.71 


29 


scintillator 


on 


1.2 X 10^ 


7 


74 - 221 


3.32 


2.01 


31 




off 


1.0 X 10^ 


5 


8-30 


0.729 


0.577 


30 




on 


1.0 X 10^ 


5 


8-30 


0.149 


0.021 


32 



Table 2: Table of some of the run parameters for the ^Li liquid scintillator and the B-plastic scintillator. The first column states whether 
the scintillator had a lead shield; the second column indicates the activity of the -^^^Cf source; the third column gives the number of runs 
that were done at different positions; the fourth column gives the range of source-detector distances r; the fifth column gives the measured 
signal rate (at the closest position only); the sixth column gives the background rate (at the closest position only); and the seventh column 
list the low-energy threshold that was used for both the data and simulation. 



The room return contribution, however, can be a signifi- 
cant fraction of the measured rate and must be taken into 
account. The room return will be very dependent upon the 
geometry and material composition of the room. In the 
CNIF the rate is largely constant over the range of mea- 
surement positions for a given source activity [39l [40] . To 
determine Rrr for each detector and each source that were 
used, measurements were taken with the source placed at 
several distances, r, from the detector. When appropri- 
ately corrected for deadtime, the intercept of a linear fit 
to the total rate in the detector versus 1/r^ gives the value 
of Rrr- Figure [6] shows the fit for data acquired using the 
^Li scintillator with the lead shield in place and using the 
higher- activity ^^^Cf source. With the measured values 
of the room return, the known source activity, and the 
calculated solid angle, one obtains the efficiency of the de- 
tector using Eq. [4j Table [3] shows the results for all eight 
configuration under which data were acquired. 

For the lifetimes, both the data and simulation were fit 
to a two-component exponential function, as discussed in 
The fast component was obtained from the 



Section 3.1 



fit to the simulation and held fixed in the fit to the data. 
For the B-plastic scintillator, the lifetime of the fast com- 
ponent was 0.47 /is, and for the ^Li-liquid scintillator it 
was equal to 1.6 jas. All the other parameters were allowed 
to float. This approach produced better fits to the data, 
particularly where the small amplitude of the component 
caused difficulty for the fit. The value of the long compo- 
nent for each configuration is given in Table [3| The stated 
uncertainties arise from statistical treatment of the data 
alone. 

3.3. Comparison with Simulation 

A Monte Carlo simulation was constructed using 
MCNP5 [41 to compare with the measurements described 
in Section |3.2[ Accurate benchmarking of this Monte 
Carlo model against the laboratory measurements is im- 
portant for understanding the response of the detector and 




Figure 6: Plot of the total neutron rate Rn versus The in- 

tercept of the linear fit gives Rrr, the contribution of neutrons that 
scatter from the room and are detected. These data are from the ^Li 
liquid scintillator with the lead annulus in place and using the higher- 
activity ^^-^Cf source. The error bars represent plus and minus the 
combine standard uncertainty. 



scintillator. The model of the test cell consisted of a right 
cylinder with the same dimensions as the actual cell and 
filled with the 0.15% ^Li-loaded liquid scintillator. The 
6 mm thick lead annulus surrounding the cylinder was also 
modeled to compare with the shielded measurements taken 
in the CNIF. The ^^^Cf source energy spectrum was as- 
sumed to be Maxwellian, and the stainless steel source 
encapsulation was modeled to take into account any at- 
tenuation arising from neutrons scattering in the source 
holder. No attempt was made to model the CNIF room- 
return or backgrounds as they were taken into account in 
the measurements. 

A neutron slowing down in a liquid scintillator loses its 
energy primarily through scattering with protons, and it 
is the interaction of each recoiling proton that produces 
the scintillation light. If the number of scintillation pho- 
tons produced is a linear function of proton recoil energy, 
the response of a detector can be modeled by converting 
the total energy deposition of a given neutron's slowing 
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Pb shield 


Activity 


Efficiency 


e (xlO-^) 


Lifetime (/is) 






(Bq) 


Data 


Simulation 


Data 


Simulation 


^Li liquid 


off 


1.2 X 10^ 


1.3 ±0.1 


1.4 ±0.1 


12.0 ±2.2 


10.2 ±0.3 


scintillator 


on 


1.2 X 10^ 


1.6 ±0.1 


1.7 ±0.1 


11.0 ± 1.9 


10.7 ± 0.2 




off 


1.0 X 10^ 


1.3 ± 0.1 


1.2 ± 0.1 


14.2 ±3.4 


10.5 ±0.3 




on 


1.0 X 10^ 


1.4 ±0.1 


1.4 ±0.1 


10.0 ±2.2 


10.8 ±0.4 


B-plastic 


off 


1.2 X 10^ 


3.0 ±0.2 


4.1 ±0.5 


1.66 ±0.20 


1.54 ±0.17 


scintillator 


on 


1.2 X 10^ 


4.0 ±0.2 


5.0 ±0.5 


1.80 ±0.07 


1.69 ±0.06 




off 


1.0 X 10^ 


3.6 ±0.2 


3.5 ±0.5 


1.78 ±0.08 


1.62 ±0.06 




on 


1.0 X 10^ 


3.9 ±0.1 


4.4 ±0.5 


1.84 ±0.07 


1.70 ±0.08 



Table 3: Comparison of the measured efficiency and capture lifetime from the MCNP model to the eight configurations of the ^Li liquid 
scintillator and B-plastic detectors. The uncertainty in the efficiency values is the standard deviation of all the individual efficiencies determined 
at each value of the detector-source distance r. The error bars for the data results come from fits to the data and represent plus and minus 
the combined standard uncertainty. The error bars for the simulation results represent the standard uncertainty. 



down history, a quantity easily obtainable from the stan- 
dard output tally of MCNP5. For incident gamma-rays 
the light-output function is linear over a wide range of en- 
ergies. For nearly all organic scintillators, however, the 
light-response function for heavy charged particles is non- 
linear. Thus, the light output from a neutron that deposits 
all of its energy in a single proton scatter is not the same 
for one that deposits all of its energy through multiple 
scatterings. To accurately model the detector response to 
fast neutrons, light output from each individual proton 
recoil in a given neutron history has to be summed to ac- 
curately model the scintillator response to neutrons. We 
assumed that the light output of proton recoils in both the 
liquid ^Li-loaded scintillator and the B-loaded plastic scin- 
tillator was the same as that for the pseudocumene-based 
liquid scintillator NE-213 and the plastic scintillator NE- 
102, which is given by 



Eee = 0.9bEp - 8.0 [1 - exp(-0.1 E^-^)] 



(5) 



where E^ is the recoiling proton's energy, and E^e is the 
electron-equivalent light output in units of MeV [42l [43] . 

A custom post-processing routine was developed to 
parse the detailed neutron history output from MCNP5 
(called a PTRAC file) and to extract the pertinent quan- 
tities from each of the individual scattering events for a 
given neutron history. Only neutron histories in which the 
neutron scattered from at least one proton and was cap- 
tured on a ^Li or ^^B nucleus were examined. For each 
history in this subset, the electron-equivalent energy was 
calculated for each proton recoil and then summed to give 
the total light output for that particular incident neutron. 
The light output was summed from the start of the neu- 
tron until the neutron was captured, and in addition, the 
time between the first proton scatter and the capture was 
recorded and histogrammed. 

The comparisons between the measurements and the 
model results are listed in Table [3j The two parameters 
on which we primarily focused were the efficiency and the 
capture lifetime. The same energy and timing cuts applied 



in the analysis of the data were applied to the Monte Carlo 
results. Ten MCNP5 runs were performed for each config- 
uration (of source strength, shielding, and detector type), 
each with a different initial seed to the random number 
generator. The results for detector efficiencies and life- 
time were averaged across these 10 runs, and the reported 
error bars are the standard deviation of the mean of those 
values. We did not quantify the systematic effects but be- 
lieve that the two more significant effects arise as a result of 
uncertainty in the energy calibration of the measurements, 
which determines the value of the energy threshold, and 
knowledge of the appropriate light output function. 

The agreement between the simulation and measured 
data is good. For the test cell containing ^Li liquid scin- 
tillator, an average efficiency for detecting fast ^^^Cf neu- 
trons was on the order of 1 x 10~^, which agreed well with 
the MCNP5 Monte-Carlo models. The efficiency for the B- 
plastic was a few times higher, which we attribute to the 
larger cross section and higher concentration of the ^^B 
in the scintillator. The simulation indicates that the effi- 
ciencies for the shielded configurations were slightly higher 
than those of the unshielded configurations. This can be 
attributed to in-scattering of neutrons from the very dense 
lead shielding surrounding the scintillator cells. Slight 
variations between the weak and strong source efficien- 
cies can be attributed to the differences in source encap- 
sulation. The time distribution of the capture events also 
agrees well with the measured lifetimes. 

4. Pulse-Shape Discrimination 

The background rejection of the scintillator was very 
good due primarily to the delayed capture on ^Li, as dis- 
In addition to timing information. 



cussed in Section 2.1 



one can further discriminate between gamma rays (elec- 
tronic recoil scattering) and neutrons (nuclear recoil scat- 
tering or neutron capture) based on the pulse waveform 
shape information. Due to the different energy deposition 
mechanisms, the temporal probability density functions 
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for scintillation light creation times produced by gamma 
rays and neutrons are not the same [44]. Hence, one can 
further discriminate background gamma events from nu- 
clear recoil and neutron capture events [45l [46l [47] based 
on the observed pulse produced by each event. 

We describe a pulse-shape discrimination (PSD) method 
based on the "Matusita distance" [48j [49l |50] between an 
observed pulse and nuclear and electronic recoil templates 
for normalized event waveforms. We also consider a stan- 
dard prompt ratio method. For both methods, we esti- 
mated a discrimination threshold for accepting approxi- 
mately 50 % of the nuclear recoil events. 

4-1- Calibration Data 

We acquired digitized pulses in two calibration mea- 
surements: one with ^^^Cs that produced essentially only 
gamma ray waveforms and another with a 2.5-MeV neu- 
tron generator that produced neutrons with some small ad- 
mixture of gamma rays. The data acquisition system that 
captured the waveforms is the same as was used for the 
previous measurements. The purpose was to estimate the 
expected value of a normalized and background-corrected 
pulse (that sums to 1) for both the nuclear recoil and elec- 
tronic recoil events. We estimated the background level 
as a trimmed mean of all values of the pulse. The values 
that contributed to the trimmed mean range from 0.125 
quantile to 0.875 quantile of the distribution of observed 
pulse values. We registered each pulse according to the 
time when the pulses steeply rises from its baseline (back- 
ground) plateau to 0.3 times its maximum value. For each 
pulse, we determined the mean value of the pulse for a 
10 ns interval before this rise time. If this mean value ex- 
ceeded 0.003 times the pulse amplitude, we rejected the 
pulse. 




50 100 150 

Time (ns) 



Figure 7: Waveform templates for nuclear recoil and electronic recoil 
events determined by cluster analysis from calibration data from a 
2.5-MeV neutron source (contaminated by gamma rays) 

and the resulting pulse was divided by its integral value. 
Values of the trimmed mean at each relative time of inter- 
est between the 0.1 and 0.9 quantiles of the distribution 
were averaged. For the 2.5-MeV source, we estimated a nu- 
clear recoil template with the same robust signal averaging 
method described above. The estimated nuclear recoil and 
electronic recoil templates from the cluster analysis agree 
well with the corresponding robust signal averaging esti- 
mates. Moreover, the estimated nuclear recoil templates 
determined from start and stop pulses for the 2.5 MeV case 
were in very close agreement for the range of amplitudes 
that we attribute to neutron capture on ^Li. 

4-3. Discrimination statistics 

The Matusita distances between a normalized pulse of 
interest, p^, and the template pulses for the electronic 
recoil Pe and nuclear recoil events Pn are 



4-2. Template Estimation 

We estimated both the nuclear and electronic recoil tem- 
plates with a k- means cluster analysis [51 from calibration 
data produced by a 2.5-MeV generator, shown in Fig. [Tj 
Joint estimation of the nuclear and electronic recoil tem- 
plates is possible because this calibration data is contami- 
nated by gamma rays. Figure [7] shows the waveform tem- 
plates determined from the 2.5-MeV source. We deter- 
mined the templates by minimizing the squared Euclidean 
distance (L^) of the normalized pulses within each of the 
two clusters. We also estimated templates with a robust 
version of cluster analysis based on an distance metric. 
In this approach, within each cluster, the median value 
rather than the mean value is computed. The robust and 
non-robust cluster analysis methods yield similar template 
estimates. 

From the ^^^Cs gamma ray source, we determined an 
electronic recoil template by a robust signal averaging 
method. Each baseline-corrected pulse was normalized so 
that its maximum value was 1. At each time sample, the 
trimmed mean of all the processed pulses was computed. 



de=Y.( Vp^-VW) f (6) 

i 

and 

dn = Y^{ VPm{i) - VpJi) )^ (7) 

i 

where i is the time increment for the digitized pulse. The 
normalized pulses sum to 1. Negative values are set to 
before taking square roots in the above equations. Our 
primary PSD statistic is 

logi? = log^. (8) 

For comparison, we also computed a prompt ratio statistic 

fp = (9) 
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Figure 8: Empirical distribution of log R statistics. 

where Xp is the integrated pulse from t = to to, and Xt 
is the integrated pulse over all times. Here, we set to to 
be the time where the nuclear and electronic recoil pulses 
cross. 

For both discrimination statistics. Figure [S] and Fig- 
ure [9) we estimate an amplitude dependent discrimination 
threshold based on events that produce logi? values less 
than 0. We then formed a curve in (amplitude, \ogR) or 
(amplitude, fp) space. For each method, we sorted the 
corresponding curve data according to amplitude bins and 
determine the median amplitude and median discrimina- 
tion statistic within each bin. In sequence, we fit a mono- 
tonic regression model [52 and then a smoothing spline to 
each curve. The degrees of freedom of the smoothing spline 
were determined by cross validation ^ We determined a 
threshold for each particular amplitude by evaluating the 
smoothing spline model at that amplitude. 

The separation between the logi? statistics appears 
more dramatic than the separation between the fp statis- 
tics for the ^^^Cs and 2.5-MeV sources. Theoretically, we 
expect that the logR statistic conveys more information 
because it is based on a 201-bin representation of the ob- 
served pulse whereas the prompt ratio is based on a 2-bin 
representation of the observed pulse. A careful quantifica- 
tion of the relative performance of PSD algorithms based 
on these two statistics is a topic for further study. One 
could also form larger bins to smooth out noise before 
computing a \ogR statistic for any pulse as discussed in 
Ref. [SH . In future experiments, our digital acquisi- 
tion system will have a higher (10-bit or 12-bit) resolu- 
tion compared to the 8-bit resolution of the data shown in 
this study. This should facilitate refinement of our PSD 
techniques. In this work, we neglected to account for the 
energy dependence of the templates. In future work, we 
may account for this dependence. 

5. Summary and Conclusions 

A liquid scintillator doped with 0.15 % ^Li by weight was 
fabricated and made into a test cell. The process of mak- 
ing the scintillator does not require complicated chemical 



Figure 9: Empirical distribution of prompt ratio statistics. The 
width of the prompt time window is determined by where the nuclear 
and electronic recoil templates cross. 

techniques, and the data acquisition system and analysis 
are straightforward. The cost of the raw materials was 
not high. The cell was tested in a known field of fast neu- 
trons and the capture-gated detection method worked as 
expected. The data acquired for both the ^Li-liquid and 
the B-plastic scintillator were compared with Monte Carlo 
simulations yielding good agreement. The detection ef- 
ficiency is not large (on the order of 10~^), but that is 
primarily due to the small size of the test cell. Starting 
with a detector of this volume, the efficiency scales more 
rapidly than the volume, so a high efficiency detector need 
not be prohibitively large or expensive. The spectroscopic 
performance of the detector was not demonstrated here, 
but in future work, we will investigate the response of the 
detector to 2.5-MeV and 14-MeV monoenergetic neutrons. 

One sample has been used in this testing over the course 
of about one year and did not shown noticeable degra- 
dation. There is some loss of light yield due to the in- 
troduction of the LiCl, but it is not prohibitively large. 
Future work will focus on quantifying optical properties 
and increasing the concentration of the ^Li in the scintil- 
lator. In principle, higher concentrations would be better 
in terms of increasing the efficiency and decreasing the 
capture time, but one also has to consider the potential 
degradation of light transmission in the scintillator. We 
note that the variety of possible ^Li-loaded scintillators has 
not been explored. LiCl was chosen initially due to its rela- 
tively high solubility in water and for the ease with which 
LiCl can be produced from isotopically-enriched lithium 
carbonate. There may be other Li-containing compounds 
where one can form stable emulsions with high concen- 
trations of ^Li. Scintillators from other manufacturers 
should be investigated for durability and optical proper- 
ties. While this paper gives data on the characteristics of 
detecting fast neutrons, the scintillator is naturally a very 
efficient detector of thermal neutrons. Future work will 
examine those properties. 

Timing information based on the delay time between 
a neutron energy deposit and its capture on ^Li enables 
one to reject backgrounds with high probability. To im- 
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prove background rejection based on this timing informa- 
tion, we presented a PSD method based on the log-ratio 
of the Matusita distances between a normahzed waveform 
and template waveforms for nuclear recoil and electronic 
recoil events. In the future, we plan to quantify the rel- 
ative performance of this multidimensional method (Fig- 
ure [8| with respect to a simpler prompt ratio discrimina- 
tion method (Figure |9|. We presented a method to de- 
termine the width of the prompt time window based on 
when the template waveforms cross. We also presented 
a nonparametric method to determine an approximate 50 
percent nuclear recoil acceptance threshold for the prompt 
ratio and Matusita distance PSD methods. 

We thank C. Bass and C. Heimbach of NIST and E. 
J. Beise, H. Breuer and T. Langford of the University of 
Maryland for useful discussions. This work was supported 
in part by the American Civil Research and Development 
Foundation, grant no. RP2-2277. 
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